Topological pumping in superconducting wires with Majorana fermions 
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We study adiabatic pumping in normal metal-superconductor hybrid systems when the super- 
conductor is in symmetry class D and therefore exhibits a phase diagram with non-trivial regions 
associated with Majorana fermions. We show that topological quantization of the pumped charge 
occurs when the phase diagram has isolated trivial regions surrounded by topological ones. We pro- 
pose a set-up where this physics emerges, which is based on a spin-orbit coupled superconducting 
nanowire in a tilted Zeeman field. 
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Introduction. — A periodic modulation of two or more 
independent parameters in a quantum system may give 
rise to a dc flow of charge in the absence of any ap- 
plied bias voltage. This transport mechanism is known as 
charge pumping [1-3] and it is termed adiabatic when the 
pumping period T is much larger than any characteristic 
time scale of the system. If the system is coupled to leads 
the pumped charge can be expressed (ignoring electron 
interaction) in terms of the scattering matrix [1, 4]. Adi- 
abatic pumping (AP) is of geometrical nature and it is 
related to the Berry phase accumulated during the cyclic 
evolution. In his pioneering work, Thouless [3] showed 
that quantized AP can occur through one-dimensional 
insulating systems as a consequence of the topological 
properties of the Hamiltonian, strictly related to the exis- 
tence of gapless points enclosed by the adiabatic pumping 
cycle in parameter space [5]. The same topological argu- 
ment holds even for a finite system attached to leads [6] . 
The recent discovery of topological insulators has also 
lead to a topological classification of spin pumps [7] . 

In this work we show how to realize topological AP 
through a normal - superconductor hybrid system. When 
time- reversal symmetry is broken (i.e. the superconduc- 
tor is in symmetry class D [8]) the topological classifi- 
cation distinguishes two possible phases of the supercon- 
ductor [9, 10], depending on whether the system sup- 
ports or not mid-gap, Majorana, states. Since quasi- 
particle transport is strongly suppressed in a sufficiently 
long superconductor due to the superconducting gap, the 
pumped charge is obtained from the evolution in time of 
the normal {ree) and Andreev {rhe) components of the 
reflection matrix r of the hybrid system as [11] 
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Adiabaticity is achieved when T is much larger than the 
quasiparticle dwell time inside the superconductor, i.e. 
the time it takes before a quasiparticle is reflected or An- 
dreev reflected. When the superconductor is in symme- 
try class D and the normal lead supports a single mode. 



the reflection matrix reduces to a single amplitude, i.e. a 
quasiparticle can undergo either a normal reflection or an 
Andreev reflection. Unitarity of the scattering matrix re- 
quires that the corresponding amplitude is a simple phase 
factor e*^ [12, 13]. The pumped charge is then given by 

where the upper (lower) sign is for normal (Andreev) re- 
flection and n is an integer since scattering amplitudes 
must be single- valued. If the phase diagram in the pa- 
rameter space is simply connected and the pumping path 
can be continuously contracted to a point, the integer n 
must be zero. On the contrary, if the phase diagram is 
non-simply connected and the pumping path, running in 
the topological phase, encloses a non-topological region, 
no restriction is enforced on n. Therefore the adiabati- 
cally pumped charge can, in principle, be topological, in 
the sense that any continuous deformation of the pump- 
ing path in parameter space does not change the integer 
n [14]. In what follows we demonstrate that topological 
AP con be obtained in a (topological) superconductor, 
even relaxing the constraint of a single open channel in 
the lead by using a quantum point contact (QPC) to cou- 
ple the lead to the superconductor. We also show that 
our results are robust against disorder. 

The model. — The system we consider for a class-D 
superconductor is a spin-orbit coupled semiconductor 
nanowire in proximity to a s-wave superconductor and 
subjected to a Zeeman field [15]. In Refs. [16, 17] it has 
been shown that when the Zeeman field is aligned parallel 
to the wire (along, e.g., the x-direction) and the pairing 
amplitude is uniform, the phase diagram as a function of 
chemical potential ji and Zeeman field Vx can be divided 
into simply-connected regions characterized by a definite 
phase. On the other hand, different regions supporting 
the same phase touch at single points. If it were possible 
to expand such contact points into regions of finite area 
(see Fig. 1) we would have a topologically non-simply 
connected phase diagram, paving the way for topolog- 
ical quantized AP. This is possible by adding a small 
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FIG. 1. (color online) Color plot of the phase of the Andreev 
reflection amplitude r^e as a function of the longitudinal Zee- 
man field Vx and the chemical potential fj, in the supercon- 
ducting region. Inside colored regions the superconductor is 
in the topologically non-trivial phase while inside white re- 
gions The = and the superconductor is in the trivial phase. 
These results have been obtained with the following set of 
parameters: W/a — 10, a/t — 0.26, and Vz/t — 0.08, where a 
is the lattice constant. Inside the normal lead /i = 14 = so 
that only a single mode is propagating. 

component of the Zeeman field in a direction perpendic- 
ular to the wire (say the z-direction) , while allowing for 
a pairing amplitude that varies along the width of the 
nanowire [17, 18]. The latter situation is consistent with 
current experiments in which the proximity effect is in- 
duced by contact with a bulk superconductor from one 
side of the nanowire [19-21]. It is thus reasonable to 
assume that the amplitude of the superconducting pair- 
ing leaking in the nanowire decreases as we move away 
from the side in contact with the superconductor. We 
have considered different profiles of the pairing ampli- 
tude along the transverse j/-direction (with < y < W, 
W being the width of the wire) [22] . 

The Hamiltonian of the superconducting nanowire 
then reads 

with 

% = [{U-n)l + V-eT]6,j~[tl-ia{iy'ija''-iy,j(Ty)]6^,j) , 

where cj^. (ci^o-) creates (destroys) an electron at site i 
with spin cr, <5(ij) restricts sites i and j to be nearest- 
neighbors, cr'' (with b = x,y,z) are spin-1/2 Pauli ma- 
trices, i^ij — X ■ dij, and v'^j = y ■ dij, with = 
{jTi — rj)/\ri — rj\ being the unit vector connecting site 
j to site i. The parameters involved in % are the hop- 
ping energy t, the spin-orbit coupling strength a, and the 
chemical potential /i. The only non-zero components of 
the Zeeman field are 14 and 14 . 

Pumping in a single-mode lead. — We first consider the 
situation in which the normal lead coupled to the su- 
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FIG. 2. (color online) Cumulative pumped charge q(t) (in 
units of e) as a function of time t (in units of the pumping 
period T). The thick solid line refers to the leftmost elliptic 
path in Fig. 1, while the thick dashed line to the rightmost 
one. At the end of the cycle, the final pumped charge Q = 
q{T) is quantized. 

perconducting nanowire supports a single open channel. 
The topological invariant Af distinguishing topologically 
trivial from non-trivial phases is related to the reflection 
matrix r at the normal- metal/superconductor (NS) in- 
terface according to N' — det r [23]. For a single-mode 
normal lead this reduces to A/" = [fee |^ — |'"/te|^- When 
the superconductor is in the topological phase (A/" = — 1), 
normal reflection is absent, r^e = 0, and the mode is per- 
fectly Andreev-reflected by the Majorana fermion at the 
NS boundary, \rhe\ = 1 [13]. On the contrary, when the 
superconductor is in the trivial phase (Af — +1), the 
mode is fully reflected {\ree\ — 1). We have evaluated 
the reflection matrix by using a method which combines 
recursive Green's function algorithms [24] with a wave- 
function matching approach [25] . In Fig. 1 the supercon- 
ductor is in the non-trivial phase (|r/je| = 1) inside col- 
ored regions and the color coding depicts the phase of the 
Andreev reflection amplitude as a function of 14 and in 
the superconducting region. Inside the white regions the 
superconductor is in the topologically trivial phase and 
The vanishes. At the boundary of such white regions the 
gap closes, signaling a quantum phase transition between 
the topologically trivial and non-trivial phases. As a con- 
sequence, the topological sector of the phase diagram is 
non-simply connected and the elliptic paths marked by 
white solid lines in Fig. 1 can not be shrunk to a point 
without crossing the phase boundary. The phase of r^e is 
almost constant throughout the parameter space with the 
exception of stripes connecting different non-topological 
(white) regions, across which the phase winds by approx- 
imately 2tt. We thus expect that when the pumping 
path coincides with non-contractiblc elliptic paths of the 
type shown in Fig. 1 these resonant stripes contribute by 
roughly e to the pumped charge Q. 

To confirm this expectation we have computed the adi- 
abatic pumped charge along both paths shown in Fig. 1 
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FIG. 3. (color online) Andreev reflection \rhe\'^ as a function 
of energy E for several increasing values of the longitudinal 
Zeeman field, 0.3 < V^/t < 0.42, at fixed fi/t = 0.32 (see 
dashed line in Fig. 1). Different curves have been offset for 
clarity. The color coding is associated with the phase of rhe 
at zero energy as in Fig. 1. 

by using a gauge-invariant discretized form of Eq. (1). 
We modeled the evolution in time of the pumping pa- 
rameters as follows 

/x(i) = /l + <5A'C0s(27rt/r) 
V^it)^%+SV, sin (27: t/T) ' ^ ^ 

where p,/t = 0.52, Sfi/t = 0.195, V^/t = 0.43, and 
SV^/t = 0.17_for the left path, while fl/t = 0.935, 
Sn/t = 0.245, V^/t = 0.53, and 6V^/t = 0.22 for the right 
one. Experimentally this evolution could be realized by 
properly designing metal gates and coils to introduce 
time-dependent gate voltages and magnetic fields acting 
only on the superconducting section of the nanowire. 

In Fig. 2 we show our numerical results for the cumula- 
tive pumped charge q{t) defined as in Eq. (1) with the in- 
tegral running from to t. The solid (dashed) line refers 
to the leftmost (rightmost) path in Fig. 1. The total 
pumped charge in a cycle Q = q{T) is strictly quantized 
and equal to and — 2e, respectively, in agreement with 
the general picture described above. In both cases the 
leading contribution to Q during the cycle comes from the 
resonant stripes that connect neighboring phase bound- 
aries. Along the rightmost path both resonant stripes 
give rise to a counterclockwise rotation of the phase of 
rhe so that their corresponding contributions sum up to 
— 2e. Along the leftmost path the winding of the Andreev 
reflection amplitude phase across the resonant stripes is 
opposite and the total pumped charge vanishes. It should 
be clear that the value of the pumped charge is topolog- 
ical, i.e. it does not change under smooth deformations 
of the pumping path in parameter space. Indeed, it is 
not possible to deform the contour in order to avoid the 
crossing of the resonant stripes without intersecting the 
phase boundary lines. 

We have verified that similar resonant stripes are 
present also in different parameter ranges or with dif- 
ferent profiles of the pairing amplitude, eventually giving 
rise to a topological pumping. We have noticed that Q/e 
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FIG. 4. (color online) Cumulative pumped charge q{t) (in 
units of e) as a function of time t (in units of the pumping pe- 
riod T) in the presence of a QPC at a distance d = 60 a from 
the NS interface. At the end of the cycle the final pumped 
charge Q = q{T) is quantized. Results for different trajecto- 
ries in parameter space are denoted with distinct line styles. 
The following set of parameters has been adopted: W/a = 10, 
a/t = 0.15 and Vz/t — 0.1. In the inset we have sketched the 
setup with a QPC at a distance d from the NS boundary 
in order to restrict the number of propagating modes. In 
simulations with disorder, random potential fiuctuations are 
considered in a region L > d inside the normal lead before 
the NS interface. 

can assume any integer value, even or odd, depending 
on the details of the system under consideration. We 
stress that the presence of the resonant stripes is strictly 
connected to the non-trivial topology of the phase dia- 
gram, and thus to the existence of two different phases 
of the superconductor, one of them supporting Majorana 
fermions. 

To further support the importance of Majorana 
fermions in the pumping mechanism, we observe that the 
resonant stripes are associated with a narrowing of the 
Majorana-induced zero-energy resonance in the Andreev 
refiection, distinctive of the non-trivial phase [26, 27]. 
This is shown in Fig. 3 where we plot |r/,ep (propor- 
tional to the differential conductance of the NS interface) 
as a function of energy. Different curves are reported 
(and offset for clarity) for several values of 14 at fixed 
n/t = 0.32, corresponding to points along the dashed 
segment in Fig. 1. 

QPC geometry. — We now show that topological AP 
can be achieved in a more realistic situation in which 
the number of open channels in the normal electrode is 
arbitrary (so that r^e and rhe are now matrices) but a 
QPC restricts the number of modes incident on the NS 
boundary (see inset in Fig. 4). The QPC consists of a 
saddle-point constriction [13, 28] at a distance d from 
the NS interface (see Fig. 4). When the pinch-off poten- 
tial of the QPC is chosen so that only a single mode is 
transmitted, it is possible to prove (we assume the super- 
conductor to be in the topological phase) that the matrix 
rhe changes in time only through a global phase factor 
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FIG. 5. (color online) Histogram of the pumped charge 
Q/e obtained for 100 disorder realizations with L/a = 200, 
d/a = 140, and Udis/t — 0.5. Light gray (red) and dark gray 
(blue) bars refer to two different sets of parameters and dis- 
tinct pumping paths. We notice that although the value of Q 
depends on the disorder configuration, it is always quantized 
and its parity is preserved. 



while Tee remains constant, so that the pumped charge 
in a cycle is still quantized. 

This is confirmed by the numerical results in Fig. 4, 
which show the charge pumped along three different non- 
contractible paths. The pumping parameters and in 
the superconductor follow the time evolution in Eq. (3), 
while in the normal lead we fix Vx = Vx and ^ = p,. 
The three curves in Fig. 4 simply differ in the values of 
Vx, SVx, fl, and Sfj,, which identify the pumping paths 
through Eq. (3). 

The pumped charge Q is quantized to integer multiples 
of e. As in the case of a single-mode lead, the value of Q 
can be interpreted as the result of contributions from the 
crossing of resonant stripes which gives rise to a twist of 
the global phase of r/je by ~ 2tt (data not shown). Each 
time the system crosses a resonant stripe approximately 
a charge e is pumped in the normal lead, corresponding 
to the step-like features in Fig. 4. When the distance 
d is varied, the interference pattern changes in such a 
way that resonant lines could come in/out of the way of 
the pumping path. We have noticed that, although the 
value of Q/e can vary, it is always strictly quantized to 
integer values and the parity of Q/e is conserved. This 
means that single resonant lines cannot appear or dis- 
appear along the pumping cycle but they can only be 
created/destroyed in pairs or change the sign of the cor- 
responding phase winding. 

We have also investigated the effect of disorder by 
introducing a random on-site potential in a region of 
length L (> d) on the left of the NS interface (see 
Fig. 4) . Fluctuations of the on-site energy lie in the inter- 
val [— C/dis/2, C/dis/2]. We have noticed that quantization 
of the pumped charge is almost exactly preserved up to 
large values of the disorder strength (Udis/t « 5 — 6) 
provided that a single mode is transmitted through the 
QPC (see Fig. 4). As shown in Fig. 5, the details of 
a specific disorder configuration affect the value of the 
pumped charge but, given the system parameters and the 



pumping path, the parity of Q/e is insensitive to smooth 
changes of the Hamiltonian. 

Conclusions - We have shown that topological adi- 
abatic pumping occurs in a class-D superconducting 
nanowire connected to a metallic lead, provided that a 
single mode of the latter is affected by the presence of 
Majorana fermions present at the endpoints of the super- 
conducting nanowire. This is the case, for example, when 
the lead supports a single propagating mode or when 
the nanowire is coupled to the lead through a quantum 
point contact. The topological nature of pumping con- 
sists in the fact that any continuous deformation of the 
pumping path in parameter space does not change the 
charge pumped in a cycle. The necessary condition to 
achieve a finite quantized value of the pumped charge is 
that the phase diagram presents a non-simply connected 
structure, where isolated non-topological regions are sur- 
rounded by connected topological ones. This is possi- 
ble by allowing both a non-uniform pairing amplitude 
and a tilted Zeeman field. Non-contractible pumping 
paths in parameter space can thus be identified within 
the topological phase. We have furthermore verified that 
the quantization of the pumped charge is robust against 
disorder. 
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